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ABSTRACT 

We present ALMA 870-/im (345 GHz) observations of two sub-millimetre galaxies 
(SMGs) drawn from an ALMA study of the 126 sub-millimeter sources from the 
LABOCA Extended Chandra Deep Field South Survey (LESS). The ALMA data 
identify the counterparts to these previously unidentified sub-millimeter sources and 
scrcndipitously detect bright emission lines in their spectra which we show are most 
likely to be [Cn] 157.74 /im emission yielding redshifts of z = 4.42 and z = 4.44. This 
blind detection rate within the 7.5-GHz bandpass of ALMA is consistent with the 
previously derived photometric redshift distribution of SMGs and suggests a modest, 
but not dominant (<25%), tail of 870-^m selected SMGs at z >4. We find that the 
ratio of £[cii] / £fir in these SMGs is much higher than seen for similarly far- infrared- 
luminous galaxies at z ~ 0, which is attributed to the more extended gas reservoirs 
in these high-redshift ULIRGs. Indeed, in one system we show that the [Cn] emission 
shows hints of extended emission on ^ 3 kpc scales. Finally, we use the volume probed 
by our ALMA survey to show that the bright end of the [Cll] luminosity function 
evolves strongly between z = and z ~ 4.4, reflecting the increased ISM cooling in 
galaxies as a result of their higher star-formation rates. These observations demon- 
strate that even with short integrations, ALMA is able to detect the dominant fine 
structure cooling lines from high-redshift ULIRGs, measure their energetics, spatially 
resolved properties and trace their evolution with redshift. 

Key words: galaxies: starburst, galaxies: evolution, galaxies: high-redshift 
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1 INTRODUCTION 

A significant fraction of the obscured star forma tion at z ^> 
1 ari s es from the most lum inous galaxies (e.g. Bl ain et al.l 
Il999l : iLeFloc'h et alj|2009h: ultra-lum inous infrared galax- 
ies (ULIRGs: ISanders fc Mirabell Il996l ) with bolometric lu- 
minosities of ^10 12 ~ 13 Lq and star-formation rates £ 100- 
lOOOMQyr -1 . At high redshift, these galaxies are some of 
the brightest sources in the sub-millimetre wave-band and 
so are frequently called "sub-millimetre galaxies" (SMGs). 

The spectral energy distribution (SED) of the dust 
emission in these luminous dusty galaxies has a Rayleigh- 
Jeans tail that produces a negative K-correction in the sub- 
millimetre wave-band, yielding a near constant flux density- 
lu minosity dependen ce with redshift over the range z ~ 1— 
6 (|Blain et alj [20021. Although the volume density of the 
SMG population evolves rapidly with redshift, (showing a 
1000-fol d increase in the space density over the ~ 10 Gyrs to 
z ~ 2.5; IChapman et al. 2005), the volume density appears 
to subsequently decline above z ^3. This decline is not a 
result of a reduction in their apparent flux density in the 
sub-milli metre, but a real redu ction in the volume density 
of SMGs (|Wardlow et al.ll201ll ). 

However, determining the precise strength of the de- 
cline in the SMG volume density requires locating the coun- 
terparts of high-redshift SMGs. This is problematic due 
to the poor resolution of single dish sub-millimetre maps, 
which means that SMGs have to be identified through cor- 
relations between their sub-millimetre emission and that in 
other wave-bands where higher spatial resolu tion is avail- 
able ( usually the radio and / o r mid-infrared e.g.llvison et al 



unique tracers of the star-formation process in these galax- 
ies. 

The brightest and best-studied of the atomic lines in the 
far-infrared is the 2 P3/2- 2 P1/2 nne structure line of singly 
ionised carbon at 157.74-/xm (here-after [Cn]). Much of the 
[Cn] emission from galaxies arises from the warm and dense 
photo-disassociation regions (PDRs) that form on the UV- 
illuminated surface of molecular clouds, although the [Cn] 
flux from diffuse Hll regions or from c ool, diffuse interstel- 
lar gas can also be significant (e.g. iMadden et al.1 1 19931 : 
lLord et alJll996h . The [Cn] emission line therefore provides 
an indication of the gas content and the extent of the gas 
reservoir in a galaxy. Far-infrared surveys of low-redshift 
galaxies from the Kuiper Airborne Observatory and ISO 
have shown that the [Cll] line can comprise ^ 1% of the to- 
tal bolometric luminosity jStacev et al.lll99ll : lB~rauher et al.l 
l200Sl : lGracia-Carpio et alj|201lh . This bright line is thus ide- 
ally suited for deriving redshifts for obscured galaxies and 
investigating their dynamics and star-formation properties. 

Early searches for [Cll], by necessity, focused on the 
highest-redshift far-infrared sources, 2 > 4, where the [Cll] 
line is shifted int o the atmospheric windows in the sub- 
millimetre (e.g. Hvison et al.1 TI99SI ; iMaiolino et al.1 120051 . 



120091 ; IWagg et al.1 



20ld ). Most of these sources host pow- 



20021 ; IChapman et al.1 12004 llvison et all 120071 ; iPope et al 



erful AGN (as well as being ULIRGs) and it was noted that 
their [Cll] lines were weak relative to Lfir, demonstrating 
the same behavior as seen in AGN-dominated ULIRGs in 
the local Universe. However, more recent observations of 
the [Cn] emission in high redshift star-formation dominated 
ULIRGs has shown that the [Cll] emission can be as bright 
as in lo cal, low- luminosity galaxies, Lfpnl / _Lfir ~0.1- 



1% (e.g. [Hailev-Dunsheath et alJfe oOS: Ivis on et al.l l2010bl : 



2006). These identifications are probabilistic as they rely IStacev et al.l 20ld ; Valtchanov et al.l 20111 ). This has been 



on empirical correlations that have both si gnificant scat- 
ter a nd which may evolve with redshift (e.g. ICarilli fc Yunl 
2000). Moreover, the SEDs in these other wave-bands have 
positive K-corrections and hence can miss the highest red- 
shift counterparts. Indeed in sub-millimetre surveys typi- 
cally 30-50% of SMGs lack "robust" counterparts in the 
radio or mid-infrared, and these may represent an uni den- 
tifie d tail of high-redshif t SMGs (e.g. iBiggs et al.ll201ll . see 
also iLindner et al.1 l201ll ) . To circumvent this problem re- 
quires ide ntifying SMGs u s ing sub-/millimetr e interferome- 
ters (e.g. iGear et alj|200d ; iFraver et al.ll200ol ). but until re- 
cently their sensitiv ity has been too low to locate large num- 
bers of SMGs (e.g. iDannerbauer et alj|200sj ; lYounger et al.l 
120071 ; IWang et al.ll201ll ; ISmolcic et alj|2012l ). However, with 
the commissioning of the Atacama Large Millimeter Ar- 
ray (ALMA), we can now construct the large samples 
of precisely-located SMGs needed to unambiguously study 
their properties and test galaxy formation models. 

Whilst the bulk of the bolometric emission from SMGs 
is radiated through continuum emission from dust grains in 
the rest-frame far-infrared, superimposed on this are a series 
of narrow atomic and molecular emission lines. By far the 
strongest of these arise from the atomic fine-structure tran- 
sitions of Carbon, Nitrogen and Oxygen in the far-infrared, 
along with weaker, but more commonly-studied molecular 
lines visible at millimetre wavelengths (e.g. 12 CO, HCN). 
These bright emission lines are an important pathway by 
which the dense gas in these galaxies cools, and so provide a 



interpreted as due to the lower ionisation field arising from 
more widely distributed star-formation activity within these 
systems, in contrast to the compact nuclear star formatio n 
seen in low-redshift ULIRGs (e.g. ISakamoto et all [2008). 
The strength of [Cll] (and other atomic lines), and its rel- 
ative strength to the far-infrared luminosity can therefore 
be used to probe the physical properties of the interstellar 
medium in high-redshift galaxies. 

We have recently undertaken an ALMA Cycle study 
at 870/^m (345 GHz) of the 126 sub-millimeter sources lo- 
cated in the 0.5° x 0.5° LABO CA Extended Ch andra Deep 
Field South Survey ("LESS"; IWeifi et al.1 120091 ). the most 
uniform sub-millimetre survey of its kind to date. These 
ALMA data yield unambiguous identifications for a large 
fraction of the sub-millimeter sources, directly pin-pointing 
the SMG responsible for the sub-millimetre emission to 
within < 0.2" (Hodge et al. 2012, in prep), without recourse 
to statistical radio/mid-infrared associations. In this letter, 
we present ALMA observations of two of the SMGs from 
our survey for which we are able to derive their redshifts 
from serendipitous identification of the [Cll] emission line 
in the ALMA data-cubes. We use the data to measure the 
energetics of the dominant fine structure lines, search for 
spatially resolved velocity structure and to provide an esti- 
mate of the evolution of the [Cn] luminosity function with 
redshift. We adopt a cosmology with J1a =0.73, fi m =0.27, 
and -ffo = 72kms _1 Mpc -1 in which 1" corresponds to a 
physical scale of 6.7 kpc at z ~4.4. 
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Figure 1. Top: 870/^m (345 GHz) velocity integrated ALMA maps of LESS 61 and LESS 65 from the Band 7 ALMA observations in 
compact configuration (these maps include the continuum and emission lines and are naturally weighted). In both cases, we detect strong 
870/xm emission in both sub-millimeter sources associated with a single SMG. The r.m.s. noise in these maps is 0.44 and 0.42 mjy beam -1 
for LESS 61 and LESS 65 respectively. Contours denote ±3, 5, l...a. In each panel we also show the 50% contour of the synthesised 
beam which is approximately 1.8" X 1.2" in both cases. These maps are centered at a : 03 32 45.87 <5 : —28 00 23.3 and a : 03 32 52.26 
8 : -2735 26.3 for LESS 61 and LESS 65 respectively (Table 1). The inset in LESS 65 shows the velocity structure in the [Cn] in this 
source. The colour scale shows the (continuum subtracted) velocity integrated cube between —210-0 km s —1 with the blue contours at 
2, 3, 4, 5...CT. The white contours denote the velocity integrated emission 0-210kms -1 (also at 2, 3, 4, 5...<r). The spatial offset between 
the blue-shifted and red-shifted emission in these two images is 0.50 ± 0.25 " (3.3 ± 1.7kpc). If this spatial and velocity offsets represents 
rotating gas reservoir, then we estimate a dynamical mass of M^ yn ~ 3.5 X 10 10 sin 2 (i) Mq. In this inset, the solid bar shows a scale of 
1". 



2 OBSERVATIONS AND REDUCTION 

Observations of the 126 sub-millimeter sources in the LESS 
survey were obtained with ALMA at 345 GHz (Band 7) 
with a dual polarisation setup in the compact configuration, 
yielding a synthesised beam of ~ 1.8" x 1.2". Our observa- 
tions cover 7.5 GHz band-width, split between the upper- 
and lower-sidebands, 336.1-339.8 GHz and 348.1-351.9 GHz. 
The observations employed 15 antennae and data for these 
two sources were obtained between 2011 October and 2011 
November in good conditions, PWVs0.5mm. The primary 
beam of the ALMA dishes, ~ 18" FWHM at our observ- 
ing frequency, is sufficient to encompass the error-circles o f 
the SMGs from the LESS maps, £5" ijWeifi et al.ll2009l ). 
even in confused situations. The observing frequency was 
selected to match the original LABOCA discovery map to 
ensure that the ALMA and LABOCA flux densities could 
be easily compared. Each galaxy was observed for a to- 
tal of ~ 120 seconds, with phase, bandpass and flux cali- 
bration based on J0403-360, J0538-440 and Mars, respec- 
tively. The data were processed wit h the Common Astron - 
omy Software Application (CAS A; [McMullin et"afl 120071 ). 
and we constructed both velocity integrated maps and dat- 
acubes with 62.5 MHz (54.5 kms -1 ) binning. The velocity 
integrated continuum maps for the two sources discussed 
here reach noise levels of a — 0.44 and 0.42 mjy beam -1 for 
LESS 61 and LESS 65 respectively. This is a factor ~3x 



more sensitive than the original LABOCA discovery map, 
while more critically the beam is ~150 x smaller in area than 
that of LABOCA (with a corresponding reduction in the po- 
sitional uncertainty of detected sources). In the datacubes, 
the r.m.s. sensitivity is ~ 3.8 mjy beam -1 per 54.5 km s -1 
(62.5 MHz) channel. The full catalog of the ALMA SMGs 
will be published in Hodge et al. (2012, in prep), whilst the 
ALMA 870^m counts and multi-wavelength properties of 
the ALMA SMGs will be published in Karim et al. (2012, 
in prep) and Simpson et al. (2012, in prep). 



3 ANALYSIS AND DISCUSSION 

IWardlow et ail l|201ll ) derive a photometric redshift distri- 
bution for the statistically identified radio- and mid-infrared 
counterparts of the LESS sub-millimeter sources, deriving a 
median z = 2.5 ± 0.5. Taking the errors on the redshifts into 
account, we expect 3.5 ±0.5 SMGs should have far-infrared 
fine structure lines (such as [Cn] 157.74/im, [Ol] 145/im, 
[Nil] 122^m or [Om] 88 fj,m) within the 7.5 GHz bandwidth of 
ALMA. By far the brightest of these lines in local ULIRGs 
is the [Cn] emission line, with a median eq uivalent width 
^ 10 x brighter than any of the other lines (|Brauher et al.l 
|2008|) . Focusing just on [Cn], we expect 1.5±0.5 SMGs 
in our sample will have detectable [Cn] emission in the 
ALMA band-pass, with corresponding redshift ranges of 
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Table 1. Co-ordinates and Photometry 





ALESS 61.1 


ALESS 65.1 


ID a 


1033245 6 — 280025 


1(133252 4—273527 


R A ("AT, MA 1 ) 6 


PI3 32 45 87 


f!3 32 52 26 


Tipr ( AT,M h\ b 


— 9£ 00 9 <: i ^ 


— 97 9fi '3 


So™ ("T.AROCAI 


5.8 ± 1.2 


5.9 ± 1.2 


Q B7n ('AT MA\c 


A Q9 -U n 4zL 




U. OUOp.111 


< 1.2 x 10 — 4 


/QQ v 1 n — 5 


U. '-JiJVJjJiLll. 


<r Q Q y 1 0~ 5 


/qq v in -5 


0.460/im 


< 5.7 x 10 -5 


< 1.8 x 10~ 4 


n 538;;m 

^ i UUU Li 111 


< 6 9 x 10~ 5 

\J * ZJ s\ -L\J 


<r 3 6 x 10~ 4 


fl R 1 / ( TTl 
UtUJl 1 1 l.l I 


<r fi 3 x i o~ 5 

\ U. vJ -A. J.VJ 


<r 6 3 x 1 0~ 5 


0.833/im 


4.7 ±1.1 x icr 4 


<1.2x 10~ 3 


0.850/im 


8.5 ±1.0 x 10" 4 


<3.6x 10~ 4 


0.903/im 


7.3±1.6 x 10~ 4 


<2.4x 10~ 3 


2.2/mi 


3.0±0.3 x icr 3 


<2.7x 10" 3 


3.6/jm 


3.3±0.3 x 10~ 3 


2.1 ±0.2 x 10" 3 


4.5/im 


3.9±0.2 x 10" 3 


1.4 ±0.2 x 10" 3 


5.8/im 


2.7±1.0 x 10~ 3 


1.4±1.0x 10" 3 


8.0/zm 


5.4±0.7 x 10" 3 


3.0±0.8 x 10" 3 


24/im 


0.036 ±0.002 


< 0.006 


250Atm 


4.3 ±1.5 


<2.7 


350/^m 


7.4 ±1.6 


7.6 ±1.4 


500/im 


10.2 ±1.7 


10.2±1.5 


1.4 GHz 


< 0.025 


< 0.025 



Notes: All flux densities are in mjy and all limits are 3 rr. 

a LESS ID from the LABOCA catalog in IWeifi et all feOOSl) . 6 
Co-ordinates in J2000. c The ALMA flux densities have been pri- 
mary beam corrected (the primary beam corrections are a factor 
1.12 and 1.04 for ALESS 61.1 and ALESS 65.1 respectively). The 
BVRIzK+lRAC photometry comprises VIMOS (U), MUSYC 
(BVRI z), HST (zssoLp); HA WK-I (K) and IRAC 3.6-8/mi imag- 
ing (see IWardlow et aljBoTll for details). The Herschel SPIRE 
photometry is measured using archival imaging and the pho- 
tometry has been de-blended for nearby sources. The radio flux 
density limits are taken from the VLA 1.4 GHz imaging used in 
iBiees et al.l toll) . 
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Figure 2. Continuum subtracted [Cn] emission line maps of 
LESS 61.1 and LESS 65.1. These are generated by subtracting the 
continuum and then integrating the cube between ± 2 x FWHM 
of the emission line. In all panels the contours denote S / N levels 
denote ± 3, 4, 5...<r. 



z — 4.399-4.461 and z = 4.590-4.656 for the upper and lower- 
side bands. 

As an initial step in our analysis of the ALMA maps of 
the LESS sub-millimeter sources, we exploit the frequency 
coverage of our observations to search for emission lines 
in the datacubes. For each > 4-er SMG in the velocity- 
integrated ALMA maps within the primary beam, we ex- 
tract the spectra and search for emission lines by attempting 
to fit a Gaussian emission line profile to the spectrum, only 
accepting the fit if the A\ 2 provides a significant improve- 
ment (A^ 2 > 25, or ~ 5rj) over a continuum-only fit. 

3.1 Source Identification 

In two sub-millimeter sources, LESSJ033245.6-280025 
and LESSJ033252.4-273527 (hereafter ALESS 61.1 and 
ALESS 65.1, respectively following the notation in Hodge 
et al. 2012, in prep), we identify bright emission lines 
with FWHM 230 ±25 and 490 ±35 km s" 1 respectively in 
the spectra of these SMGs. The significance of the emis- 
sion lines (integrated over the line) is S / N ~ 5.3 and 7.0 
for ALESS 61.1 and ALESS 65.1 respectively. In Fig. EQwe 
show the velocity-integrated (continuum plus line) ALMA 
maps of ALESS 61.1 and ALESS 65.1 while the spectra are 
shown in Fig. [3] To highlight the emission line detections 
in both ALMA SMGs, in Fig. [2] we also show the contin- 
uum subtracted emission line maps. In this figure, the map 
is made by integrating the continuum subtracted cube over 
± 2 x FWHM of the emission line. 

In both cases, the LABOCA 870/xm source is identified 
with a single high signal-to-noise (S / N ~ 8) ALMA SMG in 
the velocity integrated maps that is located to within < 0.2" 
(after cent roiding) and 3-4 " from their nomina l LABOCA 
positions (jWeifi et al.ll2009T l. iBiggs et al l (|201lh could not 
identify robust or tentative cou nterparts for either o f these 
sub-millimeter sources (see also lWardlow et al.|[2oTll ) as no 
SMGs were detected at 1.4 GHz (3rr < 25 /xJy at 1.4 GHz). 
Both sources are very faint in the Herschel /SPIRE imaging 
at 250, 350, and 500 /jm (the SMGs are weakly detected 
after de-blending nearby sources using the 24/xm and radio 
as priors; Table 1). The ALMA maps reveal that ALESS 61.1 
is associated with a weak 24 /im source that is also visible 
in the optical-mid-infrared, but ALESS 65.1 does not have 
any counterparts in the optical or near-infrared and is only 
detected weakly in the Spitzer IRAC imaging . Together, this 
sugge st that these galaxies lie at z ^3.5 fW ardlow et al.1 
1201 if ) . With the current data, we find no evidence for strong 
gravitational lensing of either of these galaxies. 

3.2 Line Identification 

If the emission lines we identified are indeed [Cn], then 
our data indicate redshifts of 2 = 4.419 and 2 = 4.445 
for ALESS 61.1 and ALESS 65.1 respectively. The rest- 
frame equivalent widths of the emission lines are 
0.44±0.06/im and 1.12 ± 0.15 fj,m, similar to the median 
equivalent width of [Cii] in local LIRGs an d ULIRGs 
(W oJ ai] = 0.76 ± 0.06,um; iBrauher et ail l2008h . We note 
that in the 7.5 GHz spectral coverage at 345 GHz, the emis- 
sion lines contribute 10-40% of the continuum emission in 
the ve locity integrated cubes, as predicted by ISmail et al.1 
il201ll) . 
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Figure 3. ALMA spectra of ALESS 61.1 and ALESS 65.1, 
extracted at the position of the peak sub-millimetre emis- 
sion from the SMGs. From the spectra, we identify bright 
emission lines in the upper side-band which we attribute to 
[Cll] emission at 2 = 4.419 and 2 = 4.445 for ALESS 61.1 and 
ALESS 65.1, respectively. The rest-frame equivalent widths of the 
lines are 0.44±0.06/mi and 1.12 ± 0.15 fim for ALESS 61.1 and 
ALESS 65.1 respectively, similar to the median equivalent width 
of [Ciil in local LIR Gs and ULIRGs (W . [C n] = 0.76 ± 0.06 /im; 
iBrauher et al,|[2008t) . In both cases, due to their selection both 
SMGs have significant continuum emission. 



We caution that there are other possible identifications 
for these lines. In particular, the emission could also be 
[Ol]145^m (z ~4.9), [Nil] 122/xm (z ~ 6.0), [Nil] 205/xm 
(z ~ 3.1) or much higher redshift, [Olll]88/im (z ~ 8.1). 
However, the luminosities of these fine structure lines are 
expe cted to be a factor ^ 1 x fainter than the [ Ch] emission 
(e.g. IBrauher et alj|200sl ; I Walter et al"1l2009bl ; iNagao et all 
I2OI2I : iDecarli et aDl2012l ) . 

Alternatively, the emission line could also be high-J 
12 CO at lower redshift (i.e. z = 1.3, 1.6, 2.0 and 2.3 for J = 7, 
8, 9, 10 respectively). However, the implied emission line 
luminosities (L = 1-7 x 10 s Lq) are a factor ~ 10 x higher 
than the 12 CO line luminosities of any local or high-redshift 
starbursts or AGN. For example, Mrk231 has a maximum 
* 2 CO luminosity (across all 12 CO lines) of £3xl0 7 L Q 
l|van der Werf et aljliuich . whilst APM 08279, a well stud- 
ied high-redshift far- infrared luminous QSO has a 12 CO (9- 
8) luminosity Leo Six 10 s Lq. Moreover, if the line high-J 
12 CO (presumably therefore a galaxy with an AGN), then 
the gas temperature must be high (^ 300 K if the 12 CO spec- 
tral line energy distribution peaks beyond J — 7). Yet the 
characteristic dust temperature would have to be ;£ 15-25 K 
(for a dust emissivity of /3 — 1.5-2.0). In addition, if these 
SMGs are lower-redshift (z ~ 1-3) starbursts or AGN, then 



these SMG s also lie significant l y off the far-infrare d-radio 
correlation |Condon et alJll99 J : llvison et al.ll2010al ). 

Taken together, we suggest that the line is most likely a 
fine structure transition, and given the brightness and equiv- 
alent width, probably [Cll]. We note that although the red- 
shifts of these two sources are similar (An ~1500kms~ 1 ), 
they are separated by 25 arcminutes on the sky (physical dis- 
tance of 10 Mpc) and so they are not physically associated. 



3.3 Broad Band Spectral Energy Distribution 

We show in Fig. [4] the observed broad-band SED for both 
galaxies. We overlay the SEP of the well-studied z = 2.3 star- 
burst SMM J2135-0102 jSwinbank et alj|2010l ; llvison et all 
l2010bl ) redshifted to z — 4.4 and normalised to the 870/im 
photometry. This demonstrates that at z ~4.4, the 250, 
350 and 500/im colours and 1.4 GHz of these galaxies are 
consistent with the expected colours of z > 4 ULIRGs (as- 
suming the pro perties of ULIRGs at z > 4 are the same as 
those at z ~2; ISchinnerer et al.ll2008l ; ICoppin et ail 120091 ; 
iDaddi et ai]|2009l ; ICapak et al.lboo'cf r 

We fit a modified black-body to the far-infrared pho- 
tometry, deriving estimates of the far-infrared (rest-frame 
8-1000^m) luminosities of L F ir = (2.1 ± 0.4) x 10 12 L 
and Lfir = (2.0±0.4) xl0 12 L o , for ALESS 61.1 and 
ALESS 65.1 respectivel y (corresponding to star-formation 
rates of ~ 500 Mm vr" 1 ; iKennicuttll 19981 ). 

To estimate the near-infrared luminosities of these sys- 
tems, we use the rest-frame UV-mid-infrared SEDs. We fit 
elliptical, Sb, single burs t and constant star- f ormat ion rate 
spectral templates from iBruzual fc Chariot! (|2003l ) to the 
0.3-8-/xm photometry (Table 1 and Fig. [4} using hyper-Z 
jBolzonella et al.ll2000l ) at the known redshift. We allow red- 
dening of Av = 0-5 in steps of 0.2 and over-plot the best- 
fitting SED models in Fig. [4] We then calculate the rest- 
frame ff-band magnitude (observed 8 fim at z ~ 4.4), which 
is less influenced by young stars than rest-frame UV or opti- 
cal bands and is relatively unaffected by dust. We derive ab- 
solute rest-frame _ff-band magnitudes of Hab = —24.8 ± 0.2 
and Hab = -24.2 ± 0.3 for ALESS 61.1 and ALESS 65.1, re- 
spectively (Table 2). These are comparable to the average H- 
band magnitudes for radio-identified sub-m illimeter source 
coun terparts at 2 ^3 (Har = —24.1 ± 0.9; IWardlow et al.l 
2011), the bulk of which are expected to be SMGs (Smail et 
al. 2012, in prep). 

Owing to the catastrophic degeneracies between star- 
formation history, age and reddenning in SED fitting for 
dusty sources, to estimate t he stellar mass e s we adot 
a simple appr o ach. Following lHainline et ail l|2009l ) and 
IWardlow et all (|201ll ) we use the if-band magnitude to- 
gether with an average mas s-to-light ratio for a likely 
SMG star-formation history. lHainline et al.l (|2009l ) esti- 
mate a //-band mass-to-light ratio for SMGs (with burst 
and constant star- formation templates), deriving an av- 
erage Lh/M* ~3.8L©/M© (for a Salpeter IMF). This 
suggests an stellar mass for ALESS 61.1 and ALESS 65.1 
of ~ 1.5 x 10 11 and 9 x 10 10 M respectively, compa- 
rable to previous estimates for SMG s tellar masses (e.g. 



lHainline et~ai] 2009 ; Wardlow et al.ll201lh . However, we cau- 
tion that Wardlow ct al. I (|201ll ) show that the uncertainties 
in the derived spectral types and ages result in an estimated 
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Figure 4. The observed optical-radio spectral energy distribu- 
tions for ALESS 61.1 and ALESS 65.1. We adopt redshifts for 
the two galaxies of z = 4.4189 and 2 = 4.4445. We model the far- 
infrared SEDs using a single component modified black body 
dust model (dashed line) by fitting the 250-, 350-, 500-/xm lim- 
its and 870-/jm photometry. We also model the rest-frame UV- 
near-infrared photometry using HYPER- z (blue line). The red solid 
curve show the SED of the well studied z = 2.3 starburst galaxy, 
SMM J2135, at the redshift of our galaxies (normalised to the 870- 
/im flux density). This demonstrates that the far-infrared colours 
are consistent with the high-redshift nature of these galaxies. 



factor of ~ 5 x uncertainty in assumed mass-to-light ratios 
and thus stellar masses. 



3.4 The [Cn] properties 

The [Cn] luminosities of our two SMGs, Ljcn] = 1.5- 
3.2 x 10 9 Lq, are comparable to the luminosities of other 
high-redshift ULIRGs and QSOs (e.g. iMaiolino et alj|200£ ; 
Walter et ail l2009al: IWagg et al l I2OICI ; llvison et all |2010b : 
De Breuck et all l201ll ; Walter et al. 2012). We can 



crudely estimate the molecular gas mass for these galax- 
ies using the average [Cn] / 12 CO (1-0) emission line ratio 
(L[cn] I ^cofi-oi = 4400±1000) from a sample of ten 2 = 2- 
4 starbursts (e.g. Hailev-Dunsheath et al. 2010l; Ivison et all 



l2010bl ; IStacev et al.ll20ld ; iDe Breuck et al.ll201ll ). although 
we caution that there is significa nt scatter in the CO /[Cn] 
luminosity ratio (e.g. Fig. 5 of IStacev et al.l |2010| shows 
that there is a range of Lrcii] /^co(l-o) °f approxi- 
mately ldex). Nevertheless, from this compilation of high- 
redshift observations, the conversion between [Cn] line lu- 



minosity and molecular gas mass scales approximately as 
M gas = 10±2 (L C [u] I Lq), suggesting a molecular gas mass 
for these two SMGs of M gas = 1-4 x 10 10 M . This molecular 
gas mass is comparable to the media n gas mass implied fo r 
SMGs at ,z = 2, M gas = 5±lxl0 10 M Q l|Bothwell et al.ll2012T l. 
Of course, direct measurements of the low- J CO emission are 
required to provide a more reliable estimate of the gas mass 
in these systems. 

We can also combine the [Cn] and far-infrared lumi- 
nosities to derive a ratio of L[cn] / £fir = (7.1 ± 1.6) x 10 -4 
and (16.0 ± 3.7) x 10~ 4 for LESS 61.1 and LESS 65.1 respec- 
tively. As Fig. [5]shows, in local luminous star- forming galax- 
ies the ratio of L[cn] / £fir declines by a factor ~ 100 over 
three orders of magnitude in far-infrared luminosity (with 
a factor ~3x scatter at any Lfir). This "[Cn] deficit" in 
local ULIRGs has a number of possible explanations, in- 
cluding enhanced contribution from continuum emi s sion i n 
dusty, high-ionisation regions IjLuhman et al. I ll998l, 120031), 
high- ionisation effects in the dense environments l|Abel et al.l 
120091 ) or en hanced contributions to the infrared luminosity 
from AGN l|Sargsvan et alj|2012h . 

The ratio we derive for the two 2 ~ 4.4 SMGs is a fac- 
tor ~ 10 x higher than expected for 2 ~ ULIRGs with 
their far- infrared luminosity (Fig. [5}, although the high 
Lrcii] / £fir is similar to other high-r edshift galaxies of 
comparable far- infrared luminosity (e.g. IStacev et al.l |2010| : 
ICox et a.1.1 l201ll : Walter et al. 2012). This has been inter- 
preted as evidence that the molecular emission does not re- 
side in a single, compact region illuminated by an intense 
UV radiation field (as is a good approximation for local 
ULIRGs), but rather that the [Cn] reservoir is more ex- 
tended, with the high L[cn] / Lfir ratio r eflecting the lower 



density of this extended med i um (e.g. [A bel et al 



2009; 
20111 ). 



lHailev-Dunsheath et al.l |2010| ; iGracia-Carpio et alj 
This interpretation is also consistent with high-resolution 
studies of SMGs, which have shown that the dust emis- 
sion and gas reservoirs can be extended over several 



kilo-parsecs (e.g. iBiggs & Ivisonl I200S 


; lYounger et al. 


2008 


iTacconi et al. 2008; Ivison et al.ll2011; 


Swinbank et al. 


|2011) 



and studies of their mid-infrared colors and spectral proper - 
ties l|Menendez-Delmestre et"al ] |2007l; lHainline et ai1l2009l ) . 
As we will show in i]3.5l the [Cn] emission in ALESS 65.1 ap- 
pears to show velocity structure, and we determine that the 
[Cn] reservoir may be extended across ~ 3 kpc, consistent 
with this interpretation. 

Finally, we note that since the [Cn] transition is a 
primary coolant within the PDRs, it provides a probe of 
the physical co nditions of the gas and interstellar radi- 
ation field (e.g. iHollenbach fc Tielensl Il999l ; IWolfire et al.l 
2003). We can therefore use the [Cn] and far-infrared lu- 
minos ity together with the PDR models of Kaufman et al. 
i|l999l ) to place upper limits on the far-UV radiation field 
strength and characteristic density of star-forming regions 
within the ISM. The [Cii]/Lfir ratio of our two galaxies, 
L[cn] / £fir = 7.1-16 x 10~ 4 , suggests a mean interstellar 
far-UV radiation field strength (G) £ 3000x that of Milky- 



Way and H 2 density ^ IP 5 ' 5 cm 3 ( e .g.lDanielson et al.|[201ll: 



see a lso lHailev-Dunsheath et alj |2010| and Stacev et alj 



|2010| ). The far-UV field strength limit derived in this way is 
comparable to that estimated using STARBURST99, which 
suggests that a 500MQyr _1 starburst produces a lumi- 
nosity of ~ 5 x 10 45 ergs _1 at a wavelength of 1000 A or 
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Figure 5. Le/t: L[qhj /Lfir ratio as a function of the far-infrared luminosity for our tw o 2 = 4.4 ALMA SMGs compared to local star- 
form ing galaxies and ULIRGs. In the plot we also include the 2 = 4.76 LESS SMG from lPe Breuck et alj l201ll) (see also | Cop pin et al.l 
2009|). We also include a number of high-redshift starburs t s and AGN from previous studies |Malhotra et al.ll200ll; IColbert et al.lll999t 
Unger et afll200d: ISpinoglio et alj|200a ICarral et aTlll994l: iBrauher et~al]|2008l : iLuhman et al.lll998l . 120031 : IStacev et al.ll2010l: ICox et ahl 



20111) . For the local data, we calculate the median (dashed line) and scatter (grey). This figure shows that the ratio of L[cn] / ^FIR 
for high-rcdshift ULIRGs is a factor ~ 10 X higher given their far-infrared luminosities compared to those at z ~ 0. Right: The [Cn] 
luminosity function at z = 4.4 from our survey compared to 2 = 0. For the z = 4.4 luminosity function, we assume that all of the [Cn]- 
emitting galaxies in the A2 = 0.12 volume covered by our observations were detected and so we stress that these calculations yield only 
a lower-limit on the vol ume density of high-r edshift [Cn] -emitters. The 2 = observations are derived from volume density of IRAS 
sources at z < 0.05 from IBrauher et al.l J200ST) ■ The red dotted line shows the predicted [Co] local luminosity function for a constant 
£[<~itt ] / Lyib. = 0-002, whi lst the solid line shows the z = far-infrared luminosity function convolved with the Z/rrjn] / ^FIR function 
from IBrauher et al ] l|200Sh (with the grey region denoting the errors in the volume density of local ULIRGs and scatter in the transfer 
function). Finding two galaxies in this volume probed by our observations indicates a ^,1000 factor increase in the number density of 
luminous [Cn] emitters from z ~ 0—4.4, equivalent to a ~ 3-4 X [Cn] luminosity evolution between z = 0—4.4. 



Table 2. Physical Properties of the Galaxies 





ALESS 61.1 


ALESS 65.1 


Ac 


350.726 ± 0.034 /im 


349.073 ± 0.032 /im 


Z [CII] 


4.4189 ±0.0004 


4.4445 ±0.0005 


I[CXX] 


2.5±0.4 Jykms" 1 


5.4±0.7Jykms~ 1 


FWHM 


230 ±25 km s" 1 


470±35kms- 1 


L [CII] 


(1.5 ±0.3) X 10 9 L 


(3.2 ±0.4) x 1O 9 L 


£<FIR 


(2.1 ±0.4) x 10 12 L 


(2.0±0.4) x 10 12 L Q 


H A B 


-24.8 ±0.2 


-24.2 ±0.3 


Lh 


(5.9±1.0) x 10 11 L 


(3.4 ±0.8) x 10 11 L Q 


M* 


1.5x 10 11 M 


9. Ox 10 10 M 



Notes: The ALMA data have been primary beam corrected before 
calculating luminosities. The primary beam corrections are 1.12 
and 1.04 for ALESS 61.1 and ALESS 65.1 respectively). 

a flux density of ~ 10 3 ' 5 x that of the Milky Way for a 
source with physical extent of 3kpc (adopting a far-UV 
field strength between 912-1103A for the Milky Way of 
1.6 x 10" 3 ergs" 1 cm" 2 ). 

3.5 Internal structure 

Given the spatial resolution of our observations (~1.4" 
FWHM), we search for spatially resolved emission and 
velocity structure within these two SMGs. From the 
ALMA data of LESS 65, the total flux of the galaxy 



(•5870 = 4.24 ± 0.49 mjy) is marginally higher than the peak 
flux in the map (Speak = 3.60 ± 0.49 mjy beam" 1 ) indicat- 
ing that the source may be marginally resolved. Using the 
[Cn] emission in this galaxy, we examine whether there 
is any evidence for velocity structure by comparing the 
spatial distribution of the integrated [Cn] emission be- 
tween -210-0 km s" 1 and 0-210 km s" 1 (i.e. ±FWHM rela- 
tive to the systemic redshift in the continuum subtracted 
cube). By centroiding the two maps we derive a spatial 
offset of 3.3±1.7kpc between the redshifted and blue- 
shifted emission (Fig. [TJ. Although tentative, this spatial 
extent is similar to the [Cn] spatial extent of the bright 
z ~5 SMG HDF 850.1 where similar data are available 
I Walter et al.ll2012l ). and also to the high-J 12 CO molecu- 
lar emission FWHM (^4kpc) fro m a sample of eight SMGs 
at z ~ 2 (|Tacconi et al.l [20061 ). If the spatial offset in 
ALESS 65.1 represents rotating gas, then we estimate a mass 
Mdyn ~ 3.5 x 10 10 sin 2 (i), which is consistent with the im- 
plied gas and stellar mass estimates. 

3.6 The evolution of the [Cn] luminosity function 

Since our ALMA survey has followed up all of the LESS 
SMGs, we can use the detection rate to investigate the evo- 
lution of the [Cn] luminosity function. In the following we as- 
sume that all of the [Cn]-emitting galaxies in the Az — 0.12 
volume covered by our observations were detected and so we 
stress that these calculations yield only a lower-limit on the 
volume density of high-redshift [Cn]-emitters. 
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To estimate the 2 ~0 [Cll] luminosity function, we 
st art from the z ~ far-infrared luminosity function 
of ISanders et al . (2003) and initially assume a constant 
I/[cn] / Lfir = 0.002 (Fig. [5]). This provides a firm upper 
limit to the bright end of z ~ [Cll] luminosity function, 
given the observed decline in L[crrj / Lfir at high Lfir 
(Fig. [5). We ca n then use the correl ation of L[cn] /Lfir 
with Lfir from iBrauher et al.l (|2008l ) (see also Fig. [5} to 
provide a more reliable estimate. We show our best esti- 
mate of the z ~ [Cll] luminosity function in Fig. [5] includ- 
ing the expected dispersion due to the scatter in the local 
L[cu] I Lfir-Lfir relation. 

To assess the reliability of this z ~ [Cll] luminosity 
function we at t empt to independently derive this using the 
IBrauher et all (|200Sl ) study of [Cll] emission in a sample 
of 227 z < 0.05 galaxies from the ISO archive. The parent 
population of this study can be approximated by an IRAS 
100-^im selected sample at z < 0.05 with flux densities be- 
tween Sioo^m = 1-1000 Jy. There are ~ 11,000 IRAS sources 
within this redshift and flux density range and so w e must 
account for the incompleteness in the IBrauher et all (|2008l ) 
sample to derive the z ~0 [Cll] luminos ity function. We 
start by assuming the IBrauher et all (|2008l ) sample is a ran- 
dom subset of the parent population (although we note it 
is likely that apparently bright sources are over-represented 
in the sample, and fainter sources correspondingly under- 
represented). We then only consider IRAS sources with flux 
den sities > 10 Jy where t he fraction of galaxies observed in 
the IBrauher et al.l l|200Sl ) sample is Z 3% of the parent pop- 
ulation and then cal culate the fraction of galaxies in the 
IBrauher et al.l l|2008l ) sample compared to the number of 
IRAS galaxies at z < 0.05 in bins of 100-/im flux density 
and redshift. We use these sampling fractions to correct the 
apparent [Cll] luminosity function within z < 0.05 (a comov- 
ing volume of 0.035 Gpc 3 ) and show these data in Fig. [5] 

Whilst the local data from IBrauher" et all (|2008l ) com- 
prise a complex mix of observations derived from a number 
of studies, the data are broadly consistent with the bright 
end of the z ~ [Cll] luminosity function we derived above 
from the far-infrared luminosity function and I/[cn] / Lfir 
transfer function. 

To place the new ALMA-identified z = 4A SMGs on 
this plot, we estimate the volume covered by our observa- 
tions of ECDFS by considering that the brightest [Cll] emit- 
ters will correspond to the brightest far-infrared sources (i.e. 
SMGs). The LESS SMGs are located within a 0.5° x 0.5° re- 
gion and assuming our ALMA observations cover the [Cll] 
with a redshift range of Az = 0.12 at z ~4.4, we derive a 
comoving volume of 2.9 x 10 5 Mpc 3 . Finding two galaxies in 
this volume with [Cll] luminosities > 1 x 10 9 Lq indicates a 
£ 1000 factor increase in the number density of luminous 
[Cn] emitters from z ~ 0-4.4. This is equivalent to a ~ 3- 
4 x [Cll] luminosity evolution between z = 0-4.4 which is 
consistent with the evolution in the far-infrared luminosity 
function (factor ~ 3x in Lfir for a fix ed volume density 
of 4> = 10~ 5 Mpc~ 1 ; IWardlow et alj|201ll ) between z = and 
z = 3, although since the far-infrared luminosity function de- 
clines beyond 2 = 3, the [Cll] luminosity function may peak 
at lower redshift than probed by our observations here. 



4 SUMMARY 

We have undertaken an ALMA study of 126 sub-millimetre 
sources from the LABOCA 870^m survey of the ECDFS 
l|Weifi et al.ll2009l ). We focus here on two high-redshift SMGs 
that are precisely located by our high resolution ALMA 
continuum observations and whose ALMA spectra detect 
bright emission lines. We interpret these lines as the far- 
infrared atomic fine structure line [Cn] 157.74, indicating 
redshifts of z — 4.4 for both galaxies. We show that the ratio 
of L[cn] / Lfir is higher than comparably luminous galaxies 
at 2 = 0, but consistent with the ratio seen in other high- 
redshift ULIRGs. We interpret this as evidence that the 
molecular emission is extended and indeed, in one of our 
SMGs our data suggest that the [Cn] is resolved on £ 3 kpc 
scales. Given the volume probed by our observations, we 
show that the [Cn] luminosity function must evolve strongly 
across the 12 Gyrs between 2 = and z ~ 4. 

These results show that the wide spectral baseline 
coverage of ALMA provides the opportunity to measure 
blind redshifts of large samples of distant, obscured galaxies 
through the detection of fine structure emission lines, such 
as [Cn], in short exposure times. 
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